RIEMANN-HILBERT APPROACH TO THE SIX- VERTEX MODEL 



PAVEL BLEHER AND KARL LIECHTY 



Abstract. The six-vertex model, or the square ice model, with domain wall boundary 
conditions (DWBC) has been introduced and solved for finite n by Korepin and Izergin. 
The solution is based on the Yang-Baxter equations and it represents the free energy in 
terms of an n x n Hankel determinant. Paul Zinn- Justin observed that the Izergin-Korepin 
formula can be re-expressed in terms of the partition function of a random matrix model with 
a nonpolynomial interaction. We use this observation to obtain the large n asymptotics of 
the six-vertex model with DWBC. The solution is based on the Riemann-Hilbert approach. 
In this paper we review asymptotic results obtained in different regions of the phase diagram. 



1. Six- vertex model 

The six-vertex model, or the model of two-dimensional ice, is stated on a square lattice 
with arrows on edges. The arrows obey the rule that at every vertex there are two arrows 



(1) 



(2) 



(3) 



(4) 



(5) 



(6) 



Figure 1. The six arrow configurations allowed at a vertex. 

pointing in and two arrows pointing out. This rule is sometimes called the ice-rule. There 
are only six possible configurations of arrows at each vertex, hence the name of the model, 
see Fig. 1. 

We will consider the domain wall boundary conditions (DWBC), in which the arrows on 
the upper and lower boundaries point into the square, and the ones on the left and right 
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boundaries point out. One possible configuration with DWBC on the 4x4 lattice is shown 
on Fig. 2. 



Figure 2. An example of 4 x 4 configuration. 

The name of the square ice comes from the two-dimensional arrangement of water molecules, 
H 2 0, with oxygen atoms at the vertices of the lattice and one hydrogen atom between each 
pair of adjacent oxygen atoms. We place an arrow in the direction from a hydrogen atom 
toward an oxygen atom if there is a bond between them. Thus, as we already noticed before, 
there are two in-bound and two out-bound arrows at each vertex. 
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Figure 3. The corresponding ice model. 

For each possible vertex state we assign a weight Wi, % = 1, . . . , 6, and define, as usual, 
the partition function, as a sum over all possible arrow configurations of the product of the 
vertex weights, 

6 

Z n = £ w(a), w(a)=l[w a{x) = l[w^\ (1.1) 

arrow configurations a xEV n i=l 
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where V n is the nxn set of vertices, o~(x) G {1, . . . , 6} is the vertex configuration of a at vertex 
x, according to Fig. 1, and Ni(a) is the number of vertices of type i in the configuration 
a. The sum is taken over all possible configurations obeying the given boundary condition. 
The Gibbs measure is defined then as 

M*) = ^. (1.2) 

Our main goal is to obtain the large n asymptotics of the partition function Z n . 

In general, the six- vertex model has six parameters: the weights Wi. However, by using 
some conservation laws we can reduce these to only two parameters. Any fixed boundary 
conditions impose some conservation laws on the six- vertex model. In the case of DWBC, 
they are 

N 1 (a) = N 2 (a), N 3 (a) = N 4 (a), N 5 (a) = N 6 (a) + n. (1.3) 
This allows us to reduce to the case 

uii = w 2 = a, w 3 = 11)4 = b, w 5 = w G = c. (1.4) 

Then by using the identity, 



,2 



Z n (a, a, b, b, c, c) = c n ' ' Z n I -, -, -, -, 1, 1 ] , (1.5) 

\c c c c J 

we can reduce to the two parameters, - and - . For details on how we make this reduction, 
see, e.g., the works [lj of Allison and Reshetikhin, [1_5] of Ferrari and Spohn, and [7] of Bleher 
and Liechty. 

2. Phase diagram of the six-vertex model 
Introduce the parameter 

A= , ° . (2.1) 

2ab V ; 

The phase diagram of the six-vertex model consists of the following three regions: the ferro- 
electric phase region, A > 1; the anti-ferroelectric phase region, A < — 1; and, the disordered 
phase region, — 1 < A < 1, see, e.g., [27]. In these three regions we parameterize the weights 
in the standard way: in the ferroelectric phase region, 

a = sinh(t-7), 6 = sinh(t + 7), c = sinh(2|7|), < (7! < t, (2.2) 

in the anti-ferroelectric phase region, 

a = sinh(7 — t), b = sinh(7 + t), c = sinh(27), \t\ < 7, (2.3) 

and in the disordered phase region, 

a = sin(7 — t), b = sin(7 + t), c = sin(27), \t\ < 7. (2.4) 

The phase diagram of the model is shown on Fig. 4. 

The phase diagram and the Bethe-Ansatz solution of the six-vertex model for periodic and 
anti-periodic boundary conditions are thoroughly discussed in the works of Lieb [23j-[26j, 
Lieb, Wu [27] . Sutherland [29J, Baxter [J], Batchelor, Baxter, O'Rourke, Yung [3]. See also 
the work of Wu, Lin [32] , in which the Pfaffian solution for the six- vertex model with periodic 
boundary conditions is obtained on the free fermion line, A = 0. 



4 PAVEL BLEHER AND KARL LIECHTY 




1 a/c 



Figure 4. The phase diagram of the model, where F, AF and D mark ferro- 
electric, antiferroelectric, and disordered phases, respectively. The circular arc 
corresponds to the so-called "free fermion" line, where A = 0, and the three 
dots correspond to 1-, 2-, and 3-enumeration of alternating sign matrices. 



3. IZERGIN-KOREPIN DETERMINANTAL FORMULA 

The six-vertex model with DWBC was introduced by Korepin in J2U], who derived an 
important recursion relation for the partition function of the model. This lead to a beautiful 
determinantal formula of Izergin and Korepin [17] , for the partition function of the six- vertex 
model with DWBC. A detailed proof of this formula and its generalizations are given in the 
paper of Izergin, Coker, and Korepin [18J. When the weights are parameterized according 
to (12. 4p . the formula of Izergin-Korepin is 



where r„ is the Hankel determinant, 



ab) n2 . . 

r n , (3.1) 



d 



+k-2, 



and 



^ = det ( , (3.2) 

Ub > Ki,k<n 



m = - h ■ (3.3) 

ab 



Observe that a, b, c have different parameterizations (|2.2p - (|2.4p in different phase regions. 
An elegant derivation of the Izergin-Korepin determinantal formula from the Yang-Baxter 
equations is given in the papers of Korepin and Zinn- Justin [21] and Kuperberg [22] . 
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One of the applications of the determinantal formula is that it implies that the partition 
function r n solves the Toda equation, 

2 d 

t nt" - r n = T n+1 T n -i, n > 1, ( ) = — , (3.4) 



cf. [28]. This was used by Korepin and Zinn- Justin [2T] to derive the free energy of the 
six-vertex model with DWBC, assuming some Ansatz on the behavior of subdominant terms 
in the large n asymptotics of the free energy. 

4. The six-vertex model with DWBC and a random matrix model 

Another application of the Izergin-Korepin determinantal formula is that r n can be ex- 
pressed in terms of a partition function of a random matrix model. The relation to the 
random matrix model was obtained and used by Zinn- Justin [33J. It can be derived as 
follows. Consider first the disordered phase region. 

4.1. Disordered phase region. For the evaluation of the Hankel determinant, it is conve- 
nient to use an integral representation of the function 

<P(t)= . , (4.1) 

Y( ' sin( 7 - t) sin( 7 + 1) v ; 

namely, to write it in the form of the Laplace transform, 

/oo 
e tx m(X)dX, (4.2) 
-oo 

where 

. . sinh £(7r- 2 7) / s 

m W = 2 a ■ (4-3) 



sinh ^7r 



Then 

d 



(4.5) 



^ = J X i e tx m(X)dX, (4.4) 
and by substituting this into the Hankel determinant, (13.21) . we obtain that 

n 

Y\[e tXl m(Xi)dXi\ det(A^ +fc - 2 ) 1 < J , fc < ri 

8=1 

/n n 
Y\[e tx 'm(Xi)dXi\ det(Af" 1 ) 1 < i , fc < n J] A*" 1 . 
i=l i=l 

Consider any permutation o G S n of variables A». From the last equation we have that 

/n n 
n[e ai ^(A J )rfA l ](-l) CT det(A.f- 1 ) 1 < l , fc < n n A ^)- ( 4 - 6 ) 
4=1 4=1 

By summing over a G S n , we obtain that 

1 r n 



4=1 
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where A(A) is the Vandermonde determinant, 

A(A) = det(A J fc - 1 ) 1 < l , fc < n = l[(\ k - A,). (U 



i<k 



Equation (14. 7p expresses r n in terms of a matrix model integral. Namely, if m(x) = e v ( x \ 
then 



1 n. 



n(n-l)/2 



7T 



where the integration is over the space of n x n Hermitian matrices. The matrix model 
integral can be solved, furthermore, in terms of orthogonal polynomials. 

Introduce monic polynomials Pk{x) = x k + . . . orthogonal on the line with respect to the 
weight 

w(x) = e tx m(x), (4.10) 

so that 

POO 

P j {x)P k (x)e tx m(x)dx = h k S nm . (4.11) 



Then it follows from (14 .7p that 

fc=0 



The orthogonal polynomials satisfy the three term recurrence relation, 

xP k (x) = P k+1 (x) + Q k P n (x) + R k P k ^(x) } (4.13) 
where R„ can be found as 



see, e.g., [30J. This gives that 



Rn = ^, (4-14) 



k 

h k = h Y[R j , (4.15) 

3=1 



where 



h = !°° e tx m(x)dx = — . sm(27) (4.16) 

7-oo sm(7 + t)sm( 7 -t) V ; 

By substituting (I4.15P into (I4.12p . we obtain that 

n— 1 

KU R n~ k - ( 4 - 17 ) 



n— 1 
k=l 
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4.2. Ferroelectric phase. In the ferroelectric phase, the parameters a, 5, and c are param- 
eterized by (12. 2p . We consider the case 7 > 0, which corresponds to the region b > a + c in 
the phase diagram. The case 7 < is similar, and a and b should be exchanged in that case. 
The function <fi is the Laplace transform of a discrete measure supported on the positive 
integers: 

= . - Sinh{2l l = 4 Y e~ 2tl sinh(2 7 0. (4.18) 

YK ' sinh(t + 7) sinh(t - 7) ^ K ' J v ' 

Then, similar to (J47TJ) , we find that 

n 2 00 n 

Tn = ^ J2 A ^) 2 II i 2e ' 2tk sinh (^)] • (4-19) 

l\,...,ln=l t=l 

This is the partition function for a discrete version of a Hermitian random matrix model, 
often called a discrete orthogonal polynomial ensemble (DOPE), and can also be solved in 
terms of orthogonal polynomials. The appropriate polynomials in this case are the monic 
polynomials P n {l) = /" + ••• with the orthogonality 

00 

Pj(l)Pk(l)w{l) = h k 5 jk , w{l) = 2e~ 2tl sinh(2 7 /) = e ~ 2tl+2 ^ - e' 211 - 2 ^ . (4.20) 

1=1 

Then it follows from (14.191) that 

n-l 

r n = 2" 2 J] h k . (4.21) 



fc=0 

4.3. Critical line between disordered and ferroelectric phase. When the parameters 
a, b, and c are such that b — a = c, (so A = 1 in (I2.ip ). the Izergin-Korepin formula is not 
directly applicable. However, we may consider a limiting case of the orthogonal polynomial 
formula (I4.2ip . On the critical line 

b a , 
= 1, (4.22) 

c c 

we fix a point, 

a a — 1 6 ot + 1 

- = — , - = — I <*>!, (4.23) 

and consider the partition function 

fa-1 a - I a + 1 a + 1 \ . . 

Z » = Z » — ' — ' — • 1 ' 1 J- (424) 

Consider the limit of ( 14.21R as 



t,7->-+0, ^ oi. (4.25) 

7 



Observe that in this limit, 



a sinh(t — 7) a — 1 6 sinh(t + 7) a + 1 

c ~ sinh(27) 2 ' c ~ sinh(27) 2 ^ ' ' 
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By (Oil . flTTjl . and fP2|) . we have 



7 I 1 ® ^ & i i 

^Tl I ) ; ; ; J- ; 1 

c c c c 



2 sinh(t — 7) sinh(t + 7) 
sinh(27) 



i 2 n-l 

n 

fc=0 



(Jfe! 



(2 • 



(4.27) 



To deal with limit ( I4.25P we need to rescale the orthogonal polynomials Pjt(i). Introduce the 
rescaled variable, 



x = 2tl- 27/, 



and the rescaled limiting weight, 



w n \x 



lim f e - M+2 ^ - e -^-^) = e"* - e" ra r 



a + 1 



> 1. 



(4.28) 



(4.29) 



t,7^+0, i->a a — 1 

7 

Consider monic orthogonal polynomials Pj(x; a) satisfying the orthogonality condition, 

Pj(x; a)P k (x; a)w Q (x)dx = h ka 5j k ■ (4.30) 

To find a relation between Pk{l) and P k (x; a), introduce the monic polynomials 

P k (x)=6 k P k (x/5), 



where 

5 = 2t- 2 7 , 

and rewrite orthogonality condition (14. lip in the form 



jk, 



1=1 



which is a Riemann sum for the integral in orthogonality condition ( 14.30ft . Therefore, 



lim P k (x) = P k (x; at), 



lim 5 2k+1 h h = h kt0t . 



Thus, if we rewrite formula (I4.27P as 



(4.31) 
(4.32) 

(4.33) 
(4.34) 



a a b b 

c c c c 



2 sinh(t — 7) sinh(t + 7) 



sinh(27)5 

we can take limit ( 14.25ft . In the limit we obtain that 



n 2 n-l rtfe+l 



n 

k=0 



S 2k+1 h k 



(4.35) 



a — 1 a — 1 a + 1 at + 1 



1,1 



a + 1 



l 2 71—1 

n 

fc=0 



^k,a 



(4.36) 
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4.4. Antiferroelectric phase. In the antiferroelectric phase, the parameters a, b, and c are 
parameterized by (12. 3p . and the function 

/ \ sinh(27) , , . , 

m = — \ y ! — , t < 7 , 4.37 

sinh (7 — t) smn(7 + t) 

is the Laplace transform of a discrete measure supported on the integers: 



Mt) = sinh(2 7 ) = 2 V e 2tl ~ 2 ^ (4 38) 

K ' sinh(7 - t) sinh(7 + t) ^ 



l=—oo 

Then 

2n 2 x 



r n = — Yl m 2 X\e 2tk ~ 2lM - (4-39) 



ni 

ll,...,l n =—00 1=1 



This is again the partition function of a DOPE, and we introduce the discrete monic poly- 
nomials P n {l) = l n + . . . via the orthogonality condition 

oo 

Pj(l)Pk(l)w(l) = h k 5 ]k , w(l) = e 2 «- 2 ^L (4.40) 

l=— oo 

Then it follows from (I4.39P that 

ra-l 

r n = 2 n2 l[h k . (4.41) 

fe=0 

4.5. Critical line between the antiferroelectric and disordered phases. When the 
parameters a,b, and c are such that a + b = c, (so A = — 1 in (12. ip ). the Izergin-Korepin 
formula is not directly applicable, and we must consider a limiting case of the orthogonal 
polynomial formula (I4.4ip . On the critical line 

a b , . , . 

- + - = 1, (4.42) 

c c 

we fix a point, 

a 1 — a b 1 + a 

- = — , - = —, -1«.<1, (4.43) 
and consider the partition function 

f 1 — a 1 — a 1 + a 1 + a \ 

«. = «.(—.—.—.— .1.1)- (4-44) 

This corresponds to taking a limit of the Izergin-Korepin formula in the antiferroelectric 
phase as t, 7 — > 0, and t/7 = a. Introduce the rescaled variable, 

x = -2tl + 2 7 /, (4.45) 

and the rescaled limiting weight, 

w a {x)= lim X ~H ( 4 - 46 ) 

where 

1 + a , d d . 

r = > . (4.47) 

1 — a 
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Consider monic orthogonal polynomials Pj(x; a) satisfying the orthogonality condition, 

Pj(x; a)P k {x; a)w a (x)dx = hk, a 5jk , (4.48) 



which can be obtained from the polynomials f)4.40p by taking the appropriate scaling limit 
as t, 7 — > 0, and t/j = a. Similar to (|4.36p . we obtain 

Z „^„(^,^,^,^, M ).(i±ii)" 2 n^. (448) 



5. Large n asymptotics of Z 



n 



The asymptotic evaluation of Z n in the different regions of the phase diagram thus reduces 
to asymptotic evaluation of different systems of orthogonal polynomials. In general, this may 
be done by formulating the orthogonal polynomials as the solution to a 2 x 2 matrix valued 
Riemann-Hilbert problem as in [16]. One may then perform the steepest descent analysis 
of Deift and Zhou [13]. In the case that the weight of orthogonality is a continuous one 
on M., this analysis was performed for weights of the form exp(— nV(x)) for a very general 
class of analytic potential functions V(x) in |12] . The analysis was adapted to the case that 
the orthogonality is with respect to a discrete measure in [2] and [10]. The steepest descent 
analysis yields the following results in the different regions of the phase diagram. 

5.1. Disordered phase. 

Theorem 5.1. (See [5].j Let the weights a, b, and c, in the six-vertex model with DWBC 
be parameterized as in \2.J$ . Then, as n — >■ oo, the partition function Z n has the asymptotic 
expansion 

Z n = Cn K F n2 (l + 0(n- £ )) , e > 0, (5.1) 

where 

irab 1 27 2 

F= k = r, 5.2 

2 7 cos (V) 12 3^-27)' 

and C > is a constant. 

This proves the conjecture of Zinn- Justin, and it gives the exact value of the exponent k. 
Let us remark, that the presence of the power-like factor n K in the asymptotic expansion 
of Z n in ( 15.11) is rather unusual from the point of view of random matrix models. Also, in 
the one-cut case the usual large n asymptotics of Z n in a non-critical random matrix model 
is the so called "topological expansion", which gives Z n as an asymptotic series in powers 
of 1/n 2 (see e.g. [T4] and [B]). In this case the asymptotic expansion goes over, in general, 
non- integer inverse powers of n (see [5]). 

It is noteworthy that, as shown in asymptotic formula (15.11) remains valid on the 
borderline between the disordered and antiferroelectric phases. In this case k — A , which 
corresponds to 7 = 0. 
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5.2. Ferroelectric phase. Bleher and Liechty [7], [8] obtained the large n asymptotics of 
Z n in the ferroelectric phase, A > 1, and also on the critical line between the ferroelectric 
and disordered phases, A = 1. In the ferroelectric phase we use parameterization (12. 2p for 
a, b and c. The large n asymptotics of Z n in the ferroelectric phase is given by the following 
theorem. 

Theorem 5.2. (See [TJ.j Let the weights a, b, and c in the six-vertex model with DWBC be 
parameterized as in h2.2\) with t > 7 > 0. For any e > 0, as n — >■ 00, 



where C = 1 — e 47 , G 



3 7-t 



CG npr, 



and F = b. 



O 



(5.3) 



On the critical line between the ferroelectric and disordered phases we use the parame- 



terization b 

Z n . 



a + 1, c = 1. The main result here is the following asymptotic formula for 



Theorem 5.3. (See 
where C > 0, 

and 



.) As n — >■ 00, 



Z n = Cn K G^F n2 [l + 0(n- 1/2 ) 



k = - , G = exp 
4 



F = b. 



^ < I V . 



(5.4) 

(5.5) 
(5.6) 



Notice that in both Theorem 15.21 and in Theorem 15.31 the limiting free energy F is the 
weight b. The ground state in this phase is unique and is achieved when there is exactly 
one c-type vertex in each row and column, and the rest of the vertices are of type b. That 
is, the diagonal consists of type 5 vertices while above the diagonal all vertices are type 3 
and below all vertices are type 4. The weight of the ground state is b n (c/b) n , and thus the 
free energy in the ferroelectric phase is completely determined by the ground state. This is 
a reflection of the fact that local fluctuations from the ground state can take place only in 
a thin neighborhood of the diagonal. The conservation laws (jl.3p forbid local fluctuations 
away from the diagonal. 

5.3. Antiferroelectric phase. The large n asymptotics in the antiferroelectric phase were 
obtained by Bleher and Liechty in [9]. They are given in the following theorem. In this 
theorem $1 and $4 are the Jacobi theta functions with elliptic nome q - 
j), and the phase u is given as 



-tt 2 /27 



^see e.g. 



CO 



7i 

2 



1 + 



7 



(5.7) 



Theorem 5.4. (See $\.) Let the weights a, b, and c in the six-vertex model with DWBC be 
parameterized as in h2.ty) . As n — )■ oo ; 



C^(noj)F n (l + 0(n- L )) 



(51 
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where C > is a constant, and 

" 27*1(0;) ' ( J 

In contrast to the disordered phase, note the lack of a power like term. In contrast to 
the ferroelectric phase, notice that the free energy depends transcendentally on the weight 
of the ground state configuration. Only in the limit as 7 — > 00, which can be regarded as 
the low temperature limit, does the the weight of the ground state become dominant. For a 
discussion of this limit, see [33] . 

6. The Riemann-Hilbert approach 

All the above asymptotic results are obtained in the Riemann-Hilbert approach, but the 
concrete asymptotic analysis of the Riemann-Hilbert problem is quite different in the different 
phase regions. Let us discuss it. 

6.1. Disordered phase region. To apply the Riemann-Hilbert approach, we introduce a 
rescaled weight as 

f nx\ 

w n (x) = w [ — ) . (6.1) 

It can be written as 
where 



7 

w n {x) = e ~ nVn{x \ (6.2) 



1 sinh ( n ( 7? — 1 ) x 1 + 

V n (x) = -Cx - - In V Y / J , C = -- (6-3) 

n sinh ( 1 

The external potential V n (x) is real analytic for any finite n, but it has logarithmic singular- 
ities on the imaginary axes, which accumulate to the origin as n — > 00. In fact, the limiting 
external potential, 

lim V n (x) = V(x) = -Cx + \x\, (6.4) 

n— >oo 

is not analytic at x — 0. The Riemann-Hilbert approach developed in [5] is based on opening 
of lenses whose boundary approaches the origin as n — > 00. This turns out to be possible 
due to the fact that the density of the equilibrium measure p n (x) for the external potential 
V n (x) diverges logarithmically at the origin as n — > 00, and as a result, the jump matrix on 
the boundary of the lenses converges to the unit matrix (for details see [5]). The calculation 
of subdominant asymptotic terms in the partition function as n — > 00 is the central difficult 
part of the work [5] , and it is done by an asymptotic analysis of the solution to the Riemann- 
Hilbert problem near the turning points and near the origin. 

6.2. Ferroelectric phase region. In the ferroelectric region, the measure of orthogonality 
is a discrete one on N. To apply the Riemann-Hilbert approach to discrete orthogonal 
polynomials, we need to rescale both the weight and the lattice that supports the measure 
so that the mesh of the lattice goes to zero as n — > 00. Introduce the rescaled lattice and 
weight 

L n = N , w n (x) = e - nxil ~V (1 - e~ 4nx ) = e~ nV ^ x) , (6.5) 
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where 

V n (x) = x (l-Q-l log (1 - e- 2n <) , < ( = \ < 1. (6.6) 

Then the orthogonality condition ( 14.201) can be written as 

E/nx\ /nx\ , . , „ , 
P] ( 2t ) Pk ( 2t ) Wn{x) = kk5jk - (6 ' 7) 

Notice that, as n — > oo, V n (x) has the limit 

lim V n (x) = x(l-C), (6.8) 

71— >00 

which would indicate that, in the large n limit, the polynomials (14.201) behave as polynomials 
orthogonal on N with a simple exponential weight. These polynomials are a special case of the 
classical Meixner polynomials, and there are exact formulae for their recurrence coefficients 
(see e.g., [H]). The monic Meixner polynomials which concern us are defined from the 
orthogonality condition 

oo 

^2Qj{t)Qh(lW = h%k, q = e 2 ^ 2t , (6.9) 

i=i 

and the normalizing constants are given exactly as 

(A;!)V +1 

Up to the constant factor, Theorem 15.21 can therefore be proven by showing that hf. and 
are asymptotically close as n — > oo. More precisely, it is shown in [7] that as k — > oo, for 

any e > 0, 

h k = h* (l + O (e-" 1 -']) . (6.11) 



Q [Kl) q 



6.3. Antiferroelectric region. In the antiferroelectric region, the orthogonal polynomials 
are with respect to a discrete weight, and we rescale the weight in (14.401) and the integer 
lattice as 

L n =(^-)z, w n (x) = e- nV(x) , V(x) = \x\-(x, ( = - < 1, (6.12) 
so that the orthogonality condition (I4.40p can be written as 

E P S (^) Pk (^) Mx) = h k 5 jk . (6.13) 

The mesh of the lattice L n is 2j/n, which places an upper constraint on the equilibrium 
measure, which is the limiting distribution of zeroes of the orthogonal polynomials. This 
upper constraint is realized. The equilibrium measure, which has density p(x), is supported 
on a single interval [a, but within that interval is an interval [a 1 , ft'} on which p(x) = I/27. 
This interval is called the saturated region, and it separates the single band of support [a, (3] 
into the two analytic bands [a, a'] and [/?', 0\. Thus in effect we have a "two-cut" situation, 
which is the source of the quasi-periodic factor $4 (nco) in Theorem 15.41 
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In principle, a problem could come from the fact that the potential V(x) is not analytic at 
the origin. However, it turns out that this point of nonanalycity is always in the saturated 
region and therefore does not present a problem in the steepest descent analysis. 

As previously noted, there is no power-like term in the asymptotic formula for Z n in the 
antiferroelectric phase. The Riemann-Hilbert approach to orthogonal polynomials generally 
gives an expansion of the normalizing constants h n in inverse powers of n. In the two-cut case, 
the coefficients in this expansion may be quasi-periodic functions of n. For the orthogonal 
polynomials (14.401) it is a tedious calculation involving the Jacobi theta functions to show 
that the term of order n~ l vanishes in the expansion of h n , which then implies the absence 
of the power-like term in Z n . 
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